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Abstract 

The perturbative stability of catenoidal soap films formed between parallel, equal radii, coaxial 
rings is studied using analytical and semi-analytical methods. Using a theorem on the nature of 
eigenvalues for a class of Sturm-Liouville operators, we show that for the given boundary con- 
ditions, azimuthally asymmetric perturbations are stable, while symmetric perturbations lead to 
an instability-a result demonstrated in Ben Amar et. al |7| using numerics and experiment. Fur- 
ther, we show how to obtain the lowest real eigenvalue of perturbations, using the semi-analytical 
Asymptotic Iteration Method (AIM). Conclusions using AIM support the analytically obtained 
result as well as the results in [S]- Finally, we compute the eigenfunctions and show, pictorially, 
how the perturbed soap film evolves in time. 
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I. INTRODUCTION 



Soap films and bubbles have, over the years, been a topic of active interest both in 
pedagogy and in research, in mathematics as well as in physics. In the eighteenth and 
nineteenth century, it was Lagrange [l| and Plateau who pioneered the study of such 
minimal surfaces which eventually led to the well-known Plateau problem in mathematics. 
Plateau, in fact was also the first to perform a series of elegant experiments with soap films 
which served as a basis for future experimental and mathematical investigations. Further 
details on the science of soap films and bubbles from a physicists' viewpoint can be found 
in the well-known book of Isenberg jsl. On the mathematical front, the book by Osserman 
^, as its title suggests, provides a survey on minimal surfaces. 

Our interest in this article is focused on one specific soap film configuration-the film formed 
between a pair of parallel, coaxial, equal radii rings (see Fig. 1). Geometrically, we know 
that the surface spanned by the film is a catenoid-a minimal (zero mean curvature) surface. 
Following original work by Plateau, an extensive analysis on this class of films was done 
in 1980 by Durand [5]. More recently, some theoretical and experimental work has been 
reported in j^. 
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FIG. 1: In the left figure, a rough sketch of a catenoidal soap film suspended between two coaxial 
circular rings of radius rg = acosh(^), showiim the parameters (the distance h between the two 
rings and the minimum radius of the film a ) [g]. In the right figure, a Mathematica 7.0 generated 
catenoid (with h = 2.0 and a = 1.0) is shown. 



Finding the shape of such a soap film is a standard problem in the calculus of variations 
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8|. The shape is obtained by rotating a catenary curve about the vertical axis, to obtain a 
catenoidal surface of revolution. Mathematically, one first writes down the surface energy 
functional V[S] for the soap film configuration, given as, 

fh 

V[S] = 27Ta + r2 dz (1.1) 

Jo 

where i^z = is the surface tension and r{z) is the radius at an axial distance z from the 
ring at 2; = (the other ring is at z = h). The Euler-Lagrange equation obtained from the 
first variation is then solved to find the surface. 

An obvious and important question to ask is - what happens if we give a small perturba- 
tion about the extremal configuration? The catenoidal shape is sustained if the configuration 
is stable. If it is unstable, it collapses to two disconnected planar discs, which is also a so- 
lution of the Euler-Lag^ange equation. To understand the stability question, which was 
discussed by Plateau and later in [5|], we expand the surface energy functional about the 
extremal configuration in the following way : 

V[S] = V[So] + 5V[So] + d'V[So] + (1.2) 

The different terms in the right hand side of the above equation are the zeroth, first and 
second variation terms. Setting the first variation term ((5y[S'o]) to zero gives the Euler- 
Lagrange equation which determines the extremal configuration. The next, second variation 
term (5^y[S'o], is crucial for us because it determines the stability of the film under pertur- 
bations. 

In jsl it was shown how the second variation leads to an associated eigenvalue problem of 
the Sturm-Liouville type. Therefore, knowing the sign of the lowest eigenvalue would deter- 
mine the stability of such soap films. If the lowest eigenvalue is negative then the catenoidal 
configuration is unstable. On the other hand, a positive lowest eigenvalue confirms its sta- 
bility. Durand [Si] analyzed in detail the stability question under azimuthally symmetric 
perturbations. Much later, in 1998, Ben Amar et. al [3] studied the stability and vibra- 
tions of catenoid-shaped smectic films, both numerically and experimentally. In our work, 
we confirm the results in Q by using exclusively analytical and semi-analytical methods. 
We also trace the time evolution of the film and pictorially demonstrate our conclusions on 
stability. 
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II. EIGENVALUE EQUATION FOR AZIMUTHALLY ASYMMETRIC PERTUR- 
BATION 



The catenoidal configuration of tlie soap film suspended between two parallel, co-axial 
circular rings is described by the following equation of a catenary curve, 

z h 

r{z) = acosh( ) (2-1) 

a 2a 

where, we have considered two rings of identical radius tq = acosh(^). The azimuthally 
asymmetric perturbed configuration is described by, 

f{z,<P)=r{z)+g{z,(P) (2.2) 

where <p is the azimuthal angle in the cylindrical co-ordinate system. g{z, (f)) is the az- 
imuthally asymmetric perturbation with boundary condition: ^^(0, (p) = g{h, 0) = and 

The surface energy functional for this case comes from the following formula p] for any 
general surface area (using cylindrical coordinates and f = {f{z, (/)) cos0, f{z, 0) sin0, z)), 



Thus we have, 

'^dcPdz (2.4) 



V[S]=2aJ^ ^ ^ ^Pif^ + l) + fl 

where, fz = ^ and f<j> = Hence, in this case, after a Taylor expansion of Eq. (12. 4p 
about the extremal configuration (Eq. (12. ip ) and performing the partial integration using 
the boundary conditions mentioned earlier, the second variation of the surface energy cor- 
responding to the extremal surface 5*0 becomes [s], 

5^V[So] = a r r ^" + ^' ~ (2.5) 
Jo J^uo cosh u 

where u{z) = ^ - A ^ iz[0,h] — > u[-uo,uo],uo = ^). The equations (I23D,(I23D and fl23|) 



are also verified in detail using the general discussions on minimal surfaces given in [9[. 
g{z, (j)) is related to the infinitesimal displacement ^{u) normal to the extremal surface of 
the film, by the following equation given in [sj , 

g{u, 0) = i{u, (j)) cosh(u) (2.6) 



Eq. (12. 6p can be rewritten using g{u, (p) = g{u, (j) + 27r) as, 



g{uA) = au) ' ^' \ cosh(ti) (2.7) 



sm mc 



where m = 0, 1, 2, 3 

The eigenvalue problem related to the second variation (Eq. (12.51) ). as given in [5]], is, 

^ — dv^ ^^"^ " cosh^(^) ^^" ^ cosh^(M)V'„ (2.8) 
where, V'n, the eigen function of Sturm-Liouville operator L, is related to C through the 
following equation, 

CX) 

C(n) = ^c„^/;„ (2.9) 

n=l 

As stated before, the stability of the soap film depends on the sign of the lowest eigenvalue 
Ai. In subsequent sections, we analyze the Eq. (12. 8p both analytically and numerically in 
order to know about the nature of the eigenvalues and the eigenf unctions. 

III. A RESULT ON STURM-LIOUVILLE OPERATORS 



We now state a theorem from the literature 10| on Sturm-Liouville operators which we use 



later to learn about the nature of the eigenvalues. 

Theorem : Let Eqs. (13.11) (see below) define a regular Sturm-Liouville problem, where 
Eq. (13. lap is the eigenvalue equation, Eqs. f l3.1bp and (I3.1cp are general boundary conditions. 



Lv := — (p(x)v')' + q{x)v = Xw{x)v, a < x < b (3.1a) 
BaV := Aiv{a) + A2v'{a) + aiv{b) + a2v'{b) = (3.1b) 
BbV := biv{a) + b2v'{a) + Biv{b) + B2v'{b) = (3.1c) 



v' = It is assumed that the interval [a, 6] is bounded (i.e. — oo < a < b < oo), the 
boundary conditions are linearly independent, all coefficients are real and the functions 
p,p',q,w are continuous in the domain [a, 6]. Also, p{x) > and w{x) > in [a,b]. If the 
boundary conditions are separated (i.e. ai = 02 = 6i = 62 = 0) and A1A2 < and B1B2 > 0, 



then, using the Green's identities 10, lUl for the Sturm-Liouville operator, we can show that 
for any eigenpair (A,f), 

^ y min[q{x) : a < x < b] 
~ 'max[w{x) : a < x < b] 
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For a proof of the above theorem see Appendix Rl 
Application: 

In our problem, x = u, p{u) = 1, q{u) = + w{u) = cosh^(M), a = —Uq, b = Uq 

and boundary conditions are separated, with A2 = B2 = and Ai,Bi 7^ 0. Therefore, 
our problem is a regular Sturm-Liouville problem satisfying the properties mentioned in the 
theorem quoted above. Thus, for Ai to be a negative eigenvalue, we need 

2 

mm[g(«) : —Uq < u < Uq] < ^ cosh^(wo) < — 5- (3.3) 

For m = 0, the above condition is always satisfied. If m = 1, we have, 1 < cosh^(Mo) < 2 . 
However, for the existence of a negative eigenvalue, Uq must be greater than 1.2, since this 
is the critical value of uq for which the lowest eigenvalue Ai is zero, for m = 0. Only above 
this value of uq does Ai become negative for m = [5]. Note that cosh^(1.2) = 3.278. Thus, 
for m = 1, a negative eigenvalue does not exist. Further, for m > 2, the inequality (13. 3p 
becomes, cosh^('Uo) < | which is impossible because cosh^(Mo) > 1 always. Therefore for all 
m > 2, no negative eigenvalues exist. 

The analysis above leads us to the result that, only the m = mode can be unstable. 



DC computed using the Asymp- 



IV. SEMI-ANALYTICAL, NUMERICAL ANALYSIS 
A. Eigenvalues using AIM 

The eigenvalues of the above Sturm-Liouville problem can 

n 

totic Iteration Method (AIM) for eigenvalue problems [121], which is briefly discussed in 
Appendix [Bl In our case the differential equation is, 

^ + (A, cosh^ u-m^ + -^^)^n = (4.1) 
du^ cosh u 

Since ipn{—U'o)= V'n(wo) = 0, ipniu) includes two factors which are {u + uq) and {u — Uq), the 

following transformation can be used for the present case, 

Ipniu) = {U^ - ul)(pn{u) (4.2) 

Using the above transformation in Eq. (14.11) . we obtain the differential equation satisfied by 

-TT = 5 2~j i^ + Kcosh u-m H r^—Wn (4.3) 

dw^ — Uq du — Uq cosh u 
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TABLE I: Comparison between the lowest eigenvalues for m = mode obtained from Asymptotic 
Iteration Method using 16 iterations in Mathematica, with the analytical values from Eqs. (j4.4p and 





Uo Uc 






Uo ^ 




Uq 


Ai (Eq. OD 


Ai (AIM) 


Uo 


Ai (Eq. 113]) 


Ai (AIM) 


1.1 


0.212 


0.244 


0.01 


24672.01 


24671.70 


1.15 


0.104 


0.107 


0.02 


6166.50 


6166.18 


1.2 


0.0 


-0.01 


0.03 


2739.56 


2739.23 


1.25 


-0.103 


-0.114 


0.04 


1540.12 


1539.80 


1.30 


-0.202 


-0.210 


0.05 


984.96 


984.64 



Comparing Eq. (14. 3 p with Eq. (jBip .using Eqs. (]B3p . f lB4p and (IBSp . we can calculate 6 (see 
Appendix B for definition) and using a particular chosen value u = Uch between — mq to uo, 
6 becomes a function of A„ and finally, A„ can be computed by finding the roots of d = 0. 
In this case, while computing the eigenvalues Uch = is used. In this way, eigenvalues were 
computed for different uq. 



B. Analytical check 

The eigenvalue equation related to this problem can be solved analytically for two limiting 
cases of the m = mode. For this azimuthally symmetric case (m = 0), when — )■ Wc = 1-2 
it is shown in that, 

Ai ^ 3.598(^^^^)i, uo<Uc (4.4a) 
2r 

Ai ^ -3.598(^^^--^)^ uo>Uc (4.4b) 
Zro 

Again if mq — )■ then, 

Ai^(^)^-2 (4.5) 
2uo 

The above check shown in Table U confirms that we can indeed use AIM to compute the 
eigenvalues for (Eq. (14.11) ). 
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TABLE II: Result of different number of iterations for finding the lowest real eigenvalue Ai for 
uq = 1.5, m = 0. 



number of iterations 



A 



9 

10 
11 
12 



-0.756028 
-0.756028 
0.448441 
0.448441 
-0.704753 



number of iterations 



A 



13 
14 
15 
16 
17 



-0.704753 
0.177329 
0.177329 
-0.660185 
-0.660185 



C. Convergence of results in AIM 

On performing the numerical calculations using AIM, we found that only 10 iterations were 
possible in Mathematica 7.0 in a 32-bit system. However, using the Improved Asymptotic 
Iteration Method (lAIM) (see Appendix[B]), 17 iterations could be carried out for the problem 
at hand. We noted that for higher value of Uq, the results using 10 iterations are different 
from the results obtained from 16 iterations. This deviation increases with increasing uq. 
The reason is that 10 iterations are not sufficient to guarantee a convergence towards the 
result. We have inspected the convergence of our result for a number of iterations choosing 
Mo = l-l- We note that the difference in results between 8th and 10th iteration is 0.06374, 
between 10th and 12th iteration it is 0.031732 , between 12th and 14th iteration it is 0.015443 
and between 14th and 16th it is 0.007601. Thus, as the iteration number increases, the result 
converges more and more. After each iteration, r/ or s (see Appendix[B]for definition) become 
zero for Uch = in an alternating fashion, resulting in the same lowest eigen value for two 
consecutive iterations. Therefore, we compare the nth and the {n + 2)th iterations. After 
16 iterations, the result is seen to converge up to the second decimal point. If we wish to 
have even better convergence, we need to further increase the number of iterations. 
However, let us inspect the same facts for uq = 1.5 (Table [IT|) . In this case the difference 
in results between 9th and 10th iteration is 1.204469, between 10th and 12th iteration is 
1.153194, between 13th, 14th iteration it is 0.882082 and between 15th and 16th iteration 
it is 0.837514. We note that the result is indeed converging as the number of iterations 
increase. But, unlike the previous case (i.e. for uq = 1.1), where after 16 iterations, the 
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result converged upto the second decimal point, we find that for uq = 1.5 the result does 
not converge even to the first decimal point! In addition, the difference between the results 
of 15th and 16th iteration is of the order of the actual eigenvalue. Clearly, for this case, we 
need a larger number of iterations to get the correct eigenvalue. 

D. A way to overcome the limitation of AIM 



In the previous section, we have seen that for mq > 1.2 we need many more iterations than 
16. However, using Mathematica we were unable to go beyond 17 iterations. In order to 
overcome this limitation, we have employed a heuristic approach which is elaborated in the 
Appendix O 

Using the combination of AIM and the above-stated method, we finally obtain the correct 
eigenvalues for m = 0. We have plotted (^)^ as a function of ^ in Fig. [2] (where , ui is the 
reduced frequency defined in 5]). Fig. |2] exactly matches with the plot shown in js] (see Fig. 
4 there). This ensures that our method is reliable and can be used to find the eigenvalues 
and analyse the stability of soap films between the two rings. 



c\j 
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FIG. 2: Plot of lowest eigenvalue(corrected) in terms of square of reduced frequency and as a 
function of 
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V. RESULT OF NUMERICAL ANALYSIS 



In this section we discuss the plots for the lowest eigenvalues as a function for different 
azimuthally asymmetric modes (m 7^ 0), where we have shown the behaviour for large Uq. 
The plots in Fig. [3] show that only the m = plot goes below the zero eigenvalue line 
(it crosses the zero eigenvalue line at Mq = 1.2 , alternatively ^ = 0.663). All other 
(m = 1,2,3,4) m 7^ curves asymptotically approach the zero value, for % — )■ 00 or 
2^ — > 0. Thus, the result of numerical analysis agrees with the fact that only the m = 
mode can be unstable, a fact which we obtained analytically as well. This is the central 
result in our article. 

5 
4 
3 

^ 2 
1 

-1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 

h/2ro 

FIG. 3: Plot of lowest eigenvalue Ai as a function of ^ for different m 

VI. PERTURBED CONFIGURATIONS 

Let us now try and see if we can obtain the perturbed configurations and visualize their 
time evolution. If an infinitesimal displacement is given to the extremal configuration of the 
soap film which is described by Eq. 12. H then the film will have a small motion about it. For 
the azimuthally symmetric case (m = 0) the general expression for the displacement about 
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;he extremal configuration, which describes this small motion of the soap film, is (given in 

3), 

oo ^ 
^{U, t) = y^ittn C0S{WJ) + — Sm{Wnt))iJniu) (6.1) 
n=l 



where 



UQ 

'ilMiu)^{u, 0) cosh^ udu (6.2) 



«0 
«0 



'0n('^)6('^) 0) cosh^ (6.3) 



Mo 

Wn = (^A„)5 = (^)^z/„ (6.4) 

and 6 = If • The frequency of oscillation Wn is related to the dimensionless reduced frequency 
Un by Eq. (16.41) . where p is the mass of the film per unit area. If Uq < ("Uq = ^ and 
Uc = 1^ = 1.2 corresponds to Ai = 0), then wl > 0, consequently wi is real and the motion 
described by Eq. (16. ip is the small oscillation of the soap film about its stable configuration. 
If uq > Uc ,wl < ,Wi is pure imaginary and the first term in Eq. ( 16. ip is hyperbolic, 

b °° 6 

^(u,t) = (ai cosh H smh.\wi\t)ipi(u) +} (a„cos(w„t)H — — sm(wnt))^niu) (6.5) 

^1 ^ Wn 

Thus, the first term grows in time and the configuration is unstable- it collapses into a 
configuration consisting of two plane discs. 

Now if we consider an azimuthally asymmetric displacement, the expression for the dis- 
placement may be written as: 

oo oo 

^{U, 0, t) = ^ ^(ttmn COs{Wmnt + m(f)) + hmn Sni{Wmnt + m(f)))i>mn{u) (6.6) 
m=0 n=l 

where 

^on = ^l I ^{u, (j),0)'ilJon{u) cosh'^ u d(p du (6.7) 

1 



uo JO 
uq ^27r 

2 



^t{u, (j), O)^on(w) cosh u d(j) du (6.8) 



1 



TT 



1 



dmn = — / / ^{u, (j),0)'ipmn{u) cosh u cos{m(f)) d(j) du (6.9) 



TT 



-uo Jo 



uo Jo 



bmn = — I I ^{u, (p,0)ijjmniu) cosh^ usm{m(j)) dcj) du (6.10) 



In Eq. ( 16. 6p . the instability occurs only in the leading order term. 
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Using the eigenvalues obtained earlier, we solve Eq. (12. 8p numerically and then plot the per- 
turbed configurations. While plotting, we treat the eigenfunction itself as the perturbation, 
though in a real situation the perturbation is the linear combination of such eigenfunctions 
satisfying the same boundary condition. 

Eq. (12. 8p . which is a second order differential equation can be written in terms of two first 
order differential equations in the following way: 



|j^ = -(Acosh'('t.)-m'^ + — fj— )C (6,12) 

du cosh (m) 

where, ( = il){u). We solve the above two ordinary differential equations numerically using 
Mathematica with the initial conditions: C,{—U{)) = and x(^'Wo)(of arbitrary choice). Using 
the computed eigenfunction we evaluate the g{u, (/), t) via the following relation: 
for the oscillating mode, 

g{u, 0, t) = ({u) cosh{u) CO s{m(j) + Wmnt) (6.13) 

for the collapsing mode, 

g{u,cl),t) = C(w) cosh('u) cosh(|woi|t) (6.14) 
Finally, we plot the time evolution of the perturbed configuration using the following relation: 

f{u,(P,t) = r{u)+g{u,(P,t) (6.15) 

In Fig. m and Fig. [5] we have shown the collapsing mode and oscillating mode respectively, 
where snapshots in time appear in each frame. The elapsed time mentioned in the figures 
are in units of ^ ^ • -^^S- El we have shown how the frequency of oscillation varies 

with the boundary (uq) and for different modes (m). 

It is clear from the time evolution that the analytical and numerical conclusions on 
stability found and stated earlier are once again confirmed. 



VII. SUMMARY 



We have shown using analytical and semi-analytical methods, that only the azimuthally 
symmetric (m = 0) perturbation of the catenoidal soap film can be unstable and this is 
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true when = ^ is greater than 1.2. Physically this means that for any given small 
m = perturbation, the soap film can collapse in an azimuthally symmetric way into two 
disconnected plane discs. On the other hand for m 7^ 0, the film remains stable as long as 
the perturbation does not become too large. This rather counter-intuitive result was shown 
first, using numerics and experiment in 7|. Our results re-confirm the result in [3] using 
different theoretical tools such clS (S-) db purely analytical method and (b) the semi-analytical 
AIM. 

One can also analyze a more general situation, where the radius of the rings are not same, 
using the methods employed here. For such a case, we have to alter the transformations in 
Eq. (14. 2 p accordingly, by invoking appropriate boundary conditions. 

In addition to the above results, our work also provides an example of the use of AIM and 
improved AIM in tackling eigenvalue problems which are difficult to solve, analytically. 
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Appendix A: Proof of the theorem referred in Section IIIII 

For a regular Sturm-Liouville problem defined by Eq. ( 13. ip . 

vLudx = —vpu'l'^^ + / {pv'u' + qvu)dx (Al) 



where, u{x),v{x) are eigenfunctions of the operator L, both satisfying the boundary condi- 
tions, Eq. (13.1 bp and Eq. (I3.1cp . Since p{x) > on [a,b], using v,the complex conjugate of 
V, in place of u in Eq. (lAip we get, 

I pb 

vLv' ^ p{a)v{a)v'{a) - p{b)v{b)v'{h) + / q\v\^ dx (A2) 



If the boundary conditions are separated, which mean, 

Aiv{a) + A2v'{a) = (A3a) 
Biv{b) + B2v'{b) = (A3b) 
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then, A1A2 ^ implies v{a)v'{a) ^ and B1B2 ^ implies v{b)v'{h) ^ 0. So, the first two 
terms on the right hand side of the inequality (IA2p are positive. Thus, we have, 

-6 pb 

vLv'dx ^ / q\v\'^dx (A4) 

J a 

For separated boundary conditions, L is self adjoint operator and hence, v and v have a 
common real eigenvalue, say A. Then, inequality ( 1A4I) becomes, 

/•b pb 

X / \v\'^wdx ^ / q\v\'^dx (A5) 



Thus, we have the theorem given by Eq. f l3.2p . 

min\q(x) : a < x < b] , , ^, 

A > ~ (A6) 

max[w[x) : a < X < b\ 

Appendix B: Asymptotic Iteration Method (AIM) 

Consider the homogeneous, linear, second order differential equation 

y" = Vo{x)y' + Soix)y (Bl) 

According to the Asymptotic Iteration Method, for sufficiently large value of n {n is an 
integer) , 

^ = ^ = a (B2) 

Vn Vn-1 

where 

Vk = v'k-i + -Sfc-i + VoVk-i (B3) 

and 

Sk = 4-1 + soVk~i (B4) 
for A; = 1, 2, 3, n. We also define a function 5 such that 

S = SnVn-l - Sn-lVn (B5) 

The eigenvalues can be computed by means of 5 = 0. 

Due to the presence of differentiation in Eqs. flB3P and flB4p . AIM may slow down the 
computer and consequently less number of iterations can be performed. However, in some 
problems a large number of iterations are required for convergence of the result. Improved 
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Asymptotic Iteration Method, introduced in [l3|, can speed up the process quite a bit. ijnix 



and Sn{x) can be expanded in series around xq, which is used in Eq. flB5l) to compute the 
eigenvalue. 

oo 

Vn{x) = ^cl{x-Xoy (B6a) 



i=0 
oo 



Sn{xo) = dl{x - xoY (B6b) 



i=0 



where and (i^'s are the Taylor's series expansion coefficients. Using Eqs. (1B6P in Eqs. (1B3P 
and flB4p we get following recursion relations (Eq. (]B7p ). 

i 

4 = (z + l)c^, + di_, + ^ c^cj;:^ (B7a) 

A;=0 

i 

d^ = {i + + doclA (B7b) 

fc=0 

<c°_, - <_ic° = (B8) 

The equation 5 = 0, now turns out to be Eq. fIBSp . which can be used for the computation 
of eigenvalues. 



Appendix C: Heuristic approach for finding eigen value 

We have used a heuristic approach to find the eigenvalues in the absence of the possibility 
of performing a large number of iterations. This approach is briefly outlined below. 

• We solve the Eq. 14.11 and plot the eigenfunction corresponding to the lowest eigenvalue, 
using the eigenvalue obtained from AIM. The method used for this has been discussed in 
Section IVll 

• We use the fact that the eigenfunction must satisfy the boundary condition i/ji^—uq) = 
'^i('^o) = and the eigenfunction will be symmetric. 

• If the plotted eigenfunction violates the above facts then we change the eigenvalue 
slightly around the value obtained from AIM and go through the same procedure until we 
get the correct eigenvalue. 

To illustrate this, let us take the example of Uq = 1.5, m = 0. The plot of the eigenfunction 
using the eigenvalue Ai = —0.66 ( which we got using AIM) and assuming {^)u=-i.5 = 

15 



0.02 has been shown in Fig. 7(a) Fig. 7(f) shows the correct eigenfunction obtained in 



the above way. The correct eigenvalue is —0.385 . 
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ai t = , initial perttiibed eeafi^ at t = 1.0 




FIG. 4: Unstable configuration ;uo = 1.5 ,m = , Ai = —0.385, x(—'Uo) = —0.1 ,unit of time 

1 fph?-]^ 
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FIG. 6: Variation of square of reduced frequency vi of oscillation with uq for different modes 




(a)Ai = -0.66 (AIM) (b)Ai = -0.8 (c)Ai = -0.5 




(d)Ai = -0.4 (e)Ai = -0.3 (f)Ai = -0.385 (correct) 



FIG. 7: Plot of ipi{u) for m = 0,uo = 1.5 taking different Ai, in order to find the correct eigenvalue. 
We start using the value that we get from AIM, which is (a), (f) is plot of the correct eigenfunction 
and therefore the correct eigenvalue is —0.385 . 
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